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Abstract 



We study the scaling behavior of the Lyapunov spectra of a chaotic shell 
model for 3D turbulence. First, we quantify localization of the Lyapunov vec- 
tors in the wavenumber space by using the numerical results. Using dimen- 
sional arguments of Kolmogorov-type, we then deduce explicitly the asymp- 
totic scaling behavior of the Lyapunov spectra. This in turn is confirmed by 
numerical results. This shell model may be regarded as a rare example of 
high-dimensional chaotic systems for which an analytic expression is known 
for the Lyapunov spectrum. Implications for the Navier-Stokes turbulence is 
given. In particular we conjecture that the distribution of Lyapunov expo- 
nents is not singular at null exponent. 
47.27.Jv, 05.45.+b 
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Typeset using REVTpjX 



We consider the behavior of Lyapunov spectra of a Gledzer's type shell model of 3D 
turbulence [1]. In this chaotic model, a set of complex variables Uj, (j — 1, 2, N), which 
represents the velocity in the shell kj = 2- 7 " 10 , are governed by the following equations of 
motion: 

j t + vkfj Uj = i[ aj u j+1 u j+2 + b j u j - 1 u j+1 + CjUj-Wj-2]* + f5 jA . (1) 

The coupling constants in the nonlinear terms are assumed as follows to ensure energy 
conservation: a,j = kj, bj = — Cj = \kj^ 2 b± = C\ = c 2 = a^-i = cln = bjy = 0. 

A number of remarkable properties were revealed by extensive numerical and analytical 
studies on this model. In particular, for large Reynolds R(= number the solution to this 
model is generally chaotic [2-3] and that its energy spectrum satisfies a Kolmogorov's scaling 
law of realistic turbulence. Other aspects of this model such as intermittency, probability 
distributions of velocity variables [4-6] and the effects of extra helicity-like invariant are also 
discussed [7,8]. 

Because these properties can be achieved with 30-50 degrees of freedom, its detailed 
study is feasible to explore possible links between the conventional theory of turbulence 
and the chaotic dynamical systems. One of the most remarkable properties is that the 
distribution function of Lyapunov exponents A appears to diverge at A = in the limit of 
large Reynolds number. This possibility was pointed out before in [9] using the f3— model. 
This suggests that the inertial subrange is connected with a large number of (almost) null 
exponents. Indeed, a correlation between Fourier and Lyapunov indexes was observed in the 
long-time average of the (squared) Lyapunov vectors. However, neither its relation to the 
Kolmogorov scaling is left unexplained nor the mechanism of accumulation of null exponents. 
In this Letter, we will see that how this characteristic Lyapunov spectrum is related with 
Kolmogorov's scaling through localization of Lyapunov vectors in the wavenumber space 
and obtain an asymptotic formula for the Lyapunov exponents. 

The equations (1) are integrated numerically by the fourth-order Runge-Kutta method, 
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together with 2N linearized equations. Numerical parameters v = 10 9 , 10 8 , 10 7 , 10 6 
respectively for N = 19, 22, 24, 27. The forcing is fixed as / = 5 x 10~ 3 x (1 + i). The time 
step used was At = 5 x 10~ 5 for N = 27. After some transient stage, the solution apparently 
reaches a statistically stationary state. Below we will consider the long-time averages over 
this state. The energy spectrum shows about 2 decades of Kolmogorov range (not shown). 
The Lyapunov dimensions are D = 19.8476, 25.4072, 30.2088, 34.9668 respectively for N = 
19,22,24,27. 

Let be the j— th Lyapunov vector for the n— th Fourier mode (j, n = 1,2, ...N). We 
plot in Fig.l the squared components of the Lyapunov vectors in time average for N = 27: 

E^\k n ) =< \v^\ 2 > . (2) 

Note that each Lyapunov vector is normalized as J2 n ,j \ v n^\ 2 = 1- Several distinct features 
are noted. 

1. Each Lyapunov vector has a support localized around a specific wavenumber. 

2. The center of the support of the Lyapunov vector lies at n rs D/2 for the largest 
Lyapunov exponent (j = 1). 

3. The central wavenumber of the support decreases with j, reaching n ~ 0. The corre- 
sponding Lyapunov exponents are positive but those corresponding to ~ are small. 

4. For larger j, the central wavenumber increases again reaching n D/2 at j D. The 
Lyapunov exponents for these are negative. 

5. For even larger j, the central wavenumber increases beyond n > D/2. These Lyapunov 
exponents asymptotically agree with the reciprocal of the viscous time scale — uk 2 of 
the equation (1). 

To summerize, for < n < D/2, there are two Lyapunov vectors for each n, one corre- 
sponding to positive exponent and the other negative. All these features are consistent with 

3 



the fact that the Lyapunov dimension measures the number of modes below the dissipation 
wavenumber and that each wavenumber has 2 degrees of freedom, the real and the imaginary 
parts of the velocity variables. 

On the basis of the above observation we will introduce the following set of hypotheses 
regarding the Lyapunov vectors in the inertial subrange j < D with D ^> 1. 

1. Each Lyapunov vector in wavenumber space is supported in a localized manner around 
a specific wavenumber. 

2. Lyapunov exponents are positive for 1 < j < D/2 and negative for D/2 < j . 

3. Let rij be the wavenumber for j— th Lyapunov vector, then 

(a) rij = D/2 - j + 1 for 1 < j < D/2, 

(b) rij = j - D/2 for D/2 < j < D. 

4. In the inertial subrange, the j— th Lyapunov exponent (j < D) is inversely proportional 
to the time scale e~ 1//3 /c~ 2 / 3 , which is characteristic to wavenumber k nj = k 2 n i . 

The last hypothesis is equivalent to assume that, through the localization of the Lyapunov 
vectors in wavenumber space, the Lyapunov exponents can be expressed solely in terms of the 
energy dissipation rate e and wavenumber k by the dimensional arguments of Kolmogorov- 
type. Combining these hypotheses, we deduce the following formulae for the Lyapunov 
exponents A^-: 



e^ 3 k 2 J 3 ~ 2- 2 ^3, forl<j<D/2 
-e 1 / 3 k 2 J 3 ~ _2 2 0'- D / 2 + 1 )/ 3 for D/2 < j £ D 



We rescale the Lyapunov exponent as \j/H, where H = Y.x^o^j is the Kolmogorov- 
Sinai entropy. This is equivalent to choose time r — at (a — YsfJa^j), such that H is 
normalized. Note that this choice of time neither influences geometric structure of strange 
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attractors nor invariant measures on them. Noting that 



we conclude that 



D/2 9 d/ 3 _ 1 



(2 2/3 - 1)^^ I 2- 2 ^ 3 (£> - oo), for l<j< D/2 

2 2(i-D/2+l)/3( 2 2/3 _ 



(5) 



for L>/2 <j~D 



22/3(2^/3 _ i) 

Because these expressions are free from arbitrary parameters, we can compare this phe- 
nomenological argument against the numerical results. 

The Lyapunov spectra were computed for four values of viscosity. In Fig. 2 we plot the 
positive Lyapunov exponents Xj/H against H together with the theoretical curve for D <C 1: 

Xj/H ~ (2 2 / 3 - 1)2- 2 ^ 3 , (6) 

which can be obtained from (5). Note that better agreement is obtained between the phe- 
nomenological theory and numerical results for larger N, that is for larger Lyapunov dimen- 
sions. In Fig. 4, we plot YjI=i Xi/H and compare with the theoretical prediction 

£Vtf =l-2-*' /s . (7) 

i=i 

Again, the theoretical prediction agree well with the numerical data and this agreement is 
better for larger N. Note that this scaling is different from the one proposed in [2] previously. 
In Fig. 4 we replotted the cumulated Lyapunov spectra using older scaling: J2i=\Xi/H vs 
j/D. We see sizable scatters in this enlarged figure and conclude that the older scaling is 
premature. 

For the other hand for j > D/2, a rough the agreement is seen between the theory and 
the numerical results. At present we do not know why the agreement is less clear. 

We discuss the implications of the above results for the realistic Navier-Stokes turbulence. 
Suppose that the support of the Lyapunov vectors are localized in the wavenumber space in 
the Navier-Stokes turbulence, as was true for the case of the shell model. This implies the 



existence of characteristic time scale in the motion represented by the Lyapunov vectors and 
that Lyapunov exponents are given by the reciprocal of relevant time scales. For simplicity, 
we consider Navier-Stokes flow under periodic boundary condition in a box of size L 3 . The 
Reynolds number is assumed to be so large that the motion we consider is in the inertial 
subrange. The total number of eddies of size r = L/n is estimated as n 3 , where n is a natural 
number. This is proportional to the number of modes whose scale is larger than r. The 
Lyapunov exponent associated with the smallest eddies can be estimated as the reciprocal 
of the Kolmogorov time-scale, that is, A ~ r ~ 2 / 3 ~ n 2 ^ 3 . In this way we find 



This shows that the number of modes whose Lyapunov exponents do not exceed A is A 9 / 2 . 
Therefore, the distribution function -P(A) of the Lyapunov exponents is 



It should be noted that this distribution -P(A) does not diverge at A = 0, in a marked 
contrast to the case of shell model. In the latter -P(A) has a singularity at A = like 1/A. In 
retrospect, The divergence of the distribution function of Lyapunov exponents stems from 
the condensation of modes at null wavenumber, which is a result of octave discretization 
of wavenumbers. Note that Ruelle's argument uses the (3— model, whose wavenumbers are 
also discretized in octaves. For realistic fluid turbulence, we conjecture that the distribution 
function of Lyapunov exponents is regular at A = because of no condensation at null 
wavenumber. Unfortunately, it is practically impossible to check the prediction (9) by 
the direct numerical simulations. Further studies on the accumulation of null Lyapunov 
exponents, e.g. in connection with the Anderson localization, will be reported elsewhere. 



n 3 ~ A 9/2 . 



(8) 



P(X) ~ A 7/2 



(9) 
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FIGURES 

FIG. 1. Time average of squared components of Lyapunov vectors : < \vn\ 2 >. Contour levels 
are 0.0489 x % % = 1, 2, 10 Two straight lines represent correspondence assumed in the hypotheses 
3. (a) and (b). 

FIG. 2. Distribution of Lyapunov exponents Xj/H. double circles (N = 19), double 
squares(iV = 22), solid squares (N = 24) and solid circles (N = 27). The dashed line denotes 
the theoretical prediction : (2 2 / 3 - 1)2~ 2 -?/ 3 . 

FIG. 3. Cumulated distribution of Lyapunov exponents J2l K/H Symbols are the same as in 
Fig.2. The dashed line denotes the theoretical prediction : 1 — 2~ 2j / 3 . 

FIG. 4. Cumulated distribution of Lyapunov exponents depicted using previous scaling: 
Y^l^i/H against j/D. Symbols are the same as in Fig.2. 

FIG. 5. Distribution of negative Lyapunov exponents \\j\/H against j — D/2. The dashed line 
denotes the theoretical prediction : (2 2 / 3 - l)/(2 2 / 3 (2 D / 3 - i))2 2 0'- D / 2 )/ 3 . 
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